Second discontinuity in the specific heat of two-phase superconductors 
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The recently proposed theoretical model of superconductors endowed with two distinct super- 
conducting phases and critical temperatures is further analyzed by introducing two distinct order 
parameters, described by two scalar fields which condensate at different temperatures. We find some 
deviations in basic thermodynamical quantities with respect to the Ginzburg-Landau one-phase su- 
perconductors. In particular, in contrast to the usual case where only one jump in specific heat 
takes place at the normal-superconductor transition temperature, we actually predict an additional 
discontinuity for Cv when passing from a superconducting phase to the other one. 

PACS numbers: 74.20.-z; 74.20.De; 11.15.Ex 

In the Ginzburg-Landau (GL) theory [3,0] the superconductivity is described in terms of a complex order parameter 
cf), which can be interpreted as the wave function of the Cooper pair in its center-of-mass frame. The classical 
phcnomcnological GL approach entails a (unique) critical temperature T@, without assuming a particular temperature- 
dependence of the coefficient a(T) appearing in the effective free energy function for unit volume, expanded up to the 
|0| 4 order: 

FcF n + a(T)|0| 2 + ^H 4 . (1) 

Quantity F n indicates the normal-phase (not superconducting) free energy density; while A, giving the strength 
of the Cooper pair binding, is assumed to be approximately constant. Ginzburg and Landau only assumed that 
the coefficient a(T) is positive above Tc, vanishes when the temperature approaches the critical value and becomes 
negative for T < Tc] around the critical temperature changes very smoothly: a(T) ~ a(Tc){T — Tc) . In the 
alternative quantum field approach analyzed below we instead adopt well-defined analytic expressions for a(T) as a 
function of the temperature. 

Actually, because of the interaction of the charged scalar field </> with the electromagnetic field A^, the order 
parameter is usually associated to the Higgs field responsible of the U(l) spontaneous symmetry breaking (SSB) 
M Si HI occurring during the normal state-superconducting phase transition. As a consequence of the symmetry 
breaking, due to a non-vanishing expectation value of the order parameter in the ground state below the critical 
temperature, the photon acquires a mass (causing the Meissncr effect) and the system becomes superconducting. By 
adopting this approach, we can initially start from a relativistically invariant Lagrangian containing the interaction 
of a single <fi with as well as the A self- interaction (hereafter h = c= 1): 

C = (D^ {D»<f>) + m^U - ±{4>U? - ± V" , (2) 

where m? > 0, F^ u = d^A v — d v A^ is the electromagnetic field strength, and = + 2ieA^ is the covariant 
derivative (2e is the electric charge of a Cooper pair). Given the above Lagrangian, the effective free energy density 
at finite temperature is formally identical to the GL expression given in |T]). However, despite the outstanding 
importance of the GL theory in superconductivity, as well as in other physical systems, it has still not been solved 
exactly beyond the mean-field approximation. Whilst this was not a serious problem for traditional superconductors, 
where the Ginzburg temperature interval is small around the critical temperature, the situation has changed especially 
with the advent of high-T c superconductors. In fact, for these systems, the Ginzburg temperature interval is large 
and we may expect strong field fluctuations and critical properties beyond the mean-field approximation. Indeed, 
in high-T c superconductors several experiments have observed critical effects in the specific heat 0, although the 
presence of a magnetic field generally makes the situation more complicated. On a theoretical side, the effect of gauge 
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field fluctuations causes great difficulties in the critical phenomena theory and, unlike the simpler <f> 4 theory for a 
neutral supcrfluid, the exact critical behaviour remains unknown. It is well-known that at the mean-field level the 
superconductive transition is discontinuous, but it seems that this result is confirmed even when field fluctuations 
are included Q . This is also confirmed by numerical simulations of lattice models for small values of the Ginzburg- 
Landau parameter x, while for large x the results are consistent with a continuous, second-order phase transition 
Q. It is thus a general belief that the standard GL model leads to a first-order transition instead of a continuous 
transition, but several other studies at one-loop (and even at two-loop) approximation have been carried out in recent 
years (see, for example, [HI U, 12, EH and references therein). Some of them entail runaway solutions of the GL 



equations (pointing towards first-order transitions), while others find a scaling behaviour with a new stable fixed 
point in the space of static parameters. Also, again beyond the plain mean-field approximation, Kleinert has shown 
the existence of a tricritical point in a superconductor, by taking the vortex fluctuations into account. In a recent 
paper Q we also considered one-loop radiative corrections to the classical GL theory and obtained, according to 
the chosen "condensation gauge" , two distinct well-defined expressions of the coefficient a as a function of T, and, 
correspondingly, two different critical temperatures T\ and T%. Let us start by expanding a complex field <f> as follows 

^i^+^e"/*, (3) 

where 770 is a real constant, 77 and 9 are real fields. Then, if we let the scalar field fluctuate around the minimum of the 
free energy, a condensation of the field 77 takes place as a result of the U(l) SSB. In Eq. the constant field t]q/\/2 
is defined as the expectation value (the condensation value) of the modulus of the scalar field cf>. Finite-temperature 
one-loop quantum corrections to the T = expression of the free energy density lead to [3] 

Fi = F a +± ai (T)rft + ±T& (4) 



with 



2 , A + 4e 2 2 ,»n 
ai = —m H — 1 . (5) 
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The parameter a\ vanishes when the temperature approaches a critical value given by 



Below T\ the expectation value of ry 2 which minimizes the free energy function results to be 



vliT) = -^. (7) 



Alternatively, we may expand the field </> as: 



*o + K + i<t>b), (8) 



where 4>o is a real constant, and (f> a , <f>b are two real scalar fields. Now we assume that a condensation takes place in 
the field (j) a (or, cquivalcntly, in <pb) rather than in the component r\. In Eq. ([8]) the constant field <j)o/y/2 is defined 
as the expectation value of the real part of (f>. In this case, after such condensation, the effective Helmholtz energy 
density writes 

Fn = F n + ±a u (T)4>l + ±<j>* (9) 

with 0] 

From the vanishing of an we now derive a different critical temperature 
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Since oo > A > 0, wc correspondingly have — Ti < Ti < T\. Accordingly, for very large self-interaction, 

A/e 2 — ► oo, we predict a maximum difference of 15% between the two critical temperatures 
Below T2 the expectation value for §\ which minimizes the free energy function is given by 

,2 f rr\ 4a n (r) / 10 N 
0o ( T ) = J — ■ ( 12 ) 

We understand the appearing of a new lower critical temperature when expanding the exponential in Eq. j3]) in 9/rjo 
and comparing with Eq. (JSj) : 

4>o ~ Vo 



,- 5 (-]+- (13) 



^ — 

The degrees of freedom carried out by the real scalar fields <j> a ,4>b are different from those corresponding to r/,0, 
and tend to coincide only in the limit tjq — ► 00. Actually, in Eqs. (|13[) the higher orders in t/q contribute at the 
denominator of the expression ([6]) as an additional A/3 term; that is an increased effective self- interaction of the 
Cooper pairs arises (A — * A ff = 4A/3) [f|. 

Since, as we have seen, two different condensations are allowed to occur inside the same system, we do not a priori 
exclude any of them. Hence we are led to introduce two order parameters, that is two scalar charged fields: the 
first one related to the condensation of the modulus of (f>\ (the corresponding phase will be hereafter denominated as 
"phase I" ) ; while the second one related to the condensation of the real part of <f>u ( "phase II" ) . 

Neglecting possible interactions between the two scalar fields, the total Lagrangian now writes: 

C = {D^ {D»<k) + m 2 4c/n - ^(4<hf + (£"Vn) + m 2 40n - ^iVn) 2 - ■ (14) 

As a matter of fact, starting from high values and then lowering the temperature we meet a first SSB at the critical 
temperature T\\ the medium becomes superconducting. Since the II-phase term au(T)(/)Q 2 + X(f>o 4 in the free energy 
density is negative for T < T2, by further lowering the temperature at T = T2 the condensation involving the second 
order-parameter is energetically favored and a new (second-order) phase transition starts. Below T2 the system is 
"more" superconducting with respect to the GL standard case since, in addition to the phase-I Cooper pairs, we 
should observe also the formation of phase-II Cooper pairs. Such two superconducting phases correspond to different 
condensations of electrons in Cooper pairs which exhibit different self-interaction, and are described by different scalar 
fields. The realization of one of the two regimes is ruled by the relative strength of the Cooper pair self-interaction 
(A) with respect to the electromagnetic interaction (e). 

Correspondingly, the total free energy density, being an additive quantity, results as the sum of contributions from 
normal-conducting electrons, phase-I superconducting Cooper pairs, and phasc-II superconducting Cooper pairs: 

F = F n for T >Ti, (15) 
F = F n + ^a 1 (T) m 2 + ^ m 4 for T 2 <T<T lt (16) 

F = F n + ^a 1 (T)r !Q 2 + ^ m i + ^an(T)cf )0 2 + ±<f )a 4 for T < T 2 , (17) 

(rjo indicates the expectation value of \<fii\; <fio indicates the expectation value of Re{0n}). 

From Eqs.Q and (|11|) wc are able to put the two free parameters of our theory, i.e. the "mass squared" m 2 and 
the self-interaction coupling constant A as functions of the two critical temperatures: 

2 T 2 T 2 

~ 2 - 6 ll12 (18) 



4(4T| - 3Tf) 



l 2e 2 (T 2 -T 2 ) 
4T 2 2 - 3T 2 
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Differences between GL and two-phases superconductors 



Therefore experimental measurements of T\ and Ti could yield an estimate of the dynamical parameters ruling the 
SSB and the electron binding in Cooper pairs. Notice that such a goal is not possible in the framework of the standard 
GL, theory where the parameters in (|TJ) are not explicitly determined. 

Inserting the above expressions in (O, 0, (fTUj) . and (fT2"]) . also the expectation values of the two scalar fields can 
be expressed in terms of T\ and Ti: 



4>o(T) 



T K T i - T 2 ) 
12(T 2 - Ti) 



(21) 



By inserting in (fT6)) and (fT7|) we may compare, for T < T2, the behavior of the free energy in the GL case, where (fl6|) 
holds also for T < T2, and in the case of two-phases superconductors for which, instead, (fT7|) applies. The free energy 
difference results to be 



AF = F, 



GL 



2ph 



e 2 T?(Tl 



T 2 ) 



2 1 2 



48(4T 2 -3T 2 )(T 2 -T 2 



(22) 



We see that such a difference increases by lowering the temperature and reaches its maximum for T = 0. 

dT 

where T = FV is the free energy. Since that the superconductive part of 

T 



dV 



The pressure is given by P 
the free energy density is independent of the volume, we have 

AP = P GL - P 2ph = -AF < . 



(23) 



Hence the pressure is expected to be larger for two-phases superconductors. Thus the differences in the free energy 
and in the pressure become more sensible far from T2, near to absolute zero. 



From S 



dF 
df 



for T < T2, we get the difference in the entropy density: 



AS = Sgl — <S2 P h = 



e 2 r 1 4 r(r 9 2 - T 2 ) 

12(4T 2 -3T 2 )(T 2 -T 2 ) 



(24) 



Being AS > 0, we can say that the two-phases superconductors are in a sense more "ordered" than the GL ones, the 
maximum difference for the entropy being reached at T = T2/V3 . 

We can compare as well the latent heat absorbed during the formation of the superconducting phase in GL and 
two-phase superconductors at a given temperature T < Ti (So indicates the entropy of the normal-phase) 



^gl(T) = T(Sq — Sgl) 
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e 2 T 2 T 2 (T 2 - T 2 ) 

12(4T 2 -3T 2 )(r 2 -T 2 ) ; (25) 



A2ph(T) — T(So — 52ph) = 

_ e 2 r 2 [T 2 4 (T 2 - T 2 ) + ^(T 2 - T 2 )] 
12(4T 2 -3T 2 )(r 2 -T 2 ) • 

The difference between Agl and A2 P h reaches its maximum at T = T2/2. 

Finally, applying the well-known formula for the specific heat at constant volume 

r -t 9S 

Cv-T — 

we obtain the difference in Cy (T < T2) 



(26) 



(27) 



Ar -r r e 2 T 1 4 T(r| - 3r 2 ) 

ACy = Cy OL - C V2ph - 12(4r 2_3 Ti2)(T 2_ r22) , (28) 

which is positive for < T < T 2 /v3, negative for T 2 j\[% < T < T 2 , and vanishes at T = T 2 j As it is seen in 
the figure, whilst AF(T 2 ) and AS(T 2 ) vanish, so that F and S are continuous in T 2 , quantity ACv{T 2 ) is not zero: 

Ag ^ 2 ) — 6( 4T 2 -^ 2 KT 2 -T 2 2 ) <0 - (29) 

We then observe a finite jump in the specific heat also in the transition from the superconducting first phase (I) to the 
second one (II), while in GL superconducting media only one discontinuity is expected (for T = T\, when the system 
changes from the normal to the superconducting regime). This sudden change in the heat capacity is a distinguishing 
characteristic of a first order phase transition. Since the jump of the specific heat in T = T\ results to be 



2 r 2 4 r 1 3 

6(4T 2 -3T 2 )(T 2 -T 2 ) 



AC V {T 1 ) = - — I, 1 2 — ^ , (30) 



the ratio between the two discontinuities can be written as 

ACy(T 2 ) = J\ 

ACV (TO T 2 



(31) 

Being the above ratio larger than 1 (and smaller than \/4/3 ), we expect the two jumps to be comparable. Thus 
also the second jump at the lower temperature — just a novel effect because it happens between two superconducting 
phases — could be experimentally investigated and measured. Notice also that, as expected, both discontinuities 
increase indefinitely in the large Cooper pairs self-interaction limit, A/e 2 — > 00, (Ti/T 2 ) 2 — > 4/3. 

A thermodynamical behavior qualitatively similar to the one predicted by our model has been recently observed 
in MgB 2 , as it is readily realized by a comparison with the existing experimental literature (see, for example (TEj 
and Refs. therein). Several two-band theories [l6| try to explain the second discontinuity in the specific heat of 
MgB 2 : probably it exists a correspondence between the two "classical-macroscopic" (since represent "collective" 
wave-functions for the condensate) order parameters in GL-likc approaches as the present one, and the two "quantum- 
microscopic" gaps in quasi-particle energy spectra predicted for MgB2 by some BCS-like theories. 

We stress that our results, as already pointed out, are a direct consequence of one-loop field theory calculations of 
the GL parameters. While computations beyond the plain mean-field approximation and their physical consequences 
have been largely considered and debated in the literature 13, 11 , 13, 13 , nevertheless those authors did not consider 



further degrees of freedom owned by a mean-field theory as the present one, which entails two concurrent order 
parameters and two differently bonded Cooper pairs. Consequently in the above-mentioned literature we do not find 
predictions of a second discontinuity in the specific heat or of some anomalous magnetic effects [13] ■ However, it is 
highly desirable that even more attention will be paid in the near future to explore the striking physical phenomena 
coming out from GL-like theories when field fluctuations are fully taken into account. 
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Finally, let us remark that attractive interactions ("Cooper-effect") and "gapped" energy spectrum characterize 
both quantum theory of fermionic superfluids (as e.g. 3 He) and BCS theory of superconductivity. Consequently we 
might expect that the basic thermodynamical properties here found for two-phase superconductors could analogously 
apply to a sort of "two-phase Fermi superfluids" (endowed with two critical temperatures) as well. 
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